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ABSTRACT 

The  asymptotic  series  expansion  for  the  motion  of  a 
charged  particle  in  a  slowly- varying  magnetic  field  has 
been  determined  by  Kruskal.   Kruskal  has  shown  that  the 
series  formally  satisfies  the  equations  of  motion.   In 
this  report  it  is  proved  that  the  series  is  in  fact  a 
valid  asymptotic  development  of  the  exact  solution. 
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NYO-7975 

ON  THE  ASYMPTOTIC  SERIES  EXPANSION  OF  THE  MOTION 
OP  A  CHARGED  PARTICLE  IN  SLOWLY- VARYING  FIELDS. 

1.   Introduction. 

The  equations  of  motion  of  a  charged  particle  In  an 
electromagnetic  field  may  be  written  as  follows: 

eR  +  By^k   -  E  =  0  (1.1) 

Here  B,  E  are  the  magnetic  and  electric  field  vectors,  and 
R  Is  the  position  vector  of  the  particle.  The  parameter  € 
Is  the  ratio 

€  =  m/e 

of  the  mass  m  of  the  particle  to  Its  charge  e.   The  dot 
signifies  differentiation  with  respect  to  the  time  t. 

Approximate  formulas  for  the  particle  motion.  I.e.  for 
solutions  of   (1.1),   were  given  by  Alfven^  '' .   Alfven's 
formulas  are  valid  If  the  field  varies  slowly.   (The  reference 
length  and  time  used  in  Alfven's  definition  of  what  is  meant 
by  a  "slow"  variation  are  the  radius  and  period  of  the  circular 
gyration  the  particle  would  execute  if  the  field  were  replaced 


(l)   H.  Alfven,  Cosmlcal  Electrodynamics,  Oxford,  The  Clarendon 
Press,  1953 
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by  a  constant  field  having  the  same  instantaneous  magnitude 
as  the  actual  field.  Thus  Alfven's  formulas  are  valid  also 
if  the  field  strength  is  large,  or  if  the  parameter  e   is 

small.)   Alfven's  results  have  been  derived  in  other  ways  and 

(2) 
have  been  generalized,  by  a  number  of  authors.   Hellwig^ 

gave  the  first  few  terms  of  an  asymptotic  series  for  the 

(^) 
particle  motion.   Kruskal^^''  has  given  the  complete 

asymptotic  series. 

We  assume  that   B  and  E  are  given  as  functions 

of  the  space  coordinates  R,   the  time  t,   and  the  parameter 

€^^^.   We  assume  that  B,   E,   and  e~  (B »  E)   are  continuous 

functions  of  R,   t,   having  continuous  derivatives  of  all 

orders  with  respect  to  R,   t.   We  assume  that  any  of  the 

components  of  B,   E,   or  €"  ( B  •  E) ,   or  any  of  their 

derivatives  with  respect  to  R,   t   is  bounded  in  absolute 


value  uniformly  in  R,   t,   €.   Also,  if   B  =\/B  •  B  is 

the  magnetic  field  strength,  we  assume  that   B~    is  bounded 

uniformly  in  R,   t,   e.   Thus   B  is  assiimed  to  have  an 

absolute  upper  bound  B  „  ,   and  also  to  have  an  absolute 
^^  max 

(2)  G.  Hellwig,  Uber  die  Bewegung  geladener  Teilchen  in  schwach 
verSnderlichen  Magnetf eldern,  Zeitschrift  filir  Naturf orschung 
7,  508-516,  1955- 

(3)  M.  D.  Kruskal,  paper  to  be  published  in  Nuovo  Cimento. 

(4)  The  notion  of  allowing   B,   E  to  vary  with   e   is 
due  to  Kruskal 
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lower  bound   B„.   >  0. 
rain 

We  are  concerned  with  the  behavior,  as   e   tends  to 
zero,  of  the  solution  of  the  equations  of  motion   (l.l) 
which  satisfies  the  initial  conditions 

for  t=0;R=R,  v=v  (1.2) 

Here  R",  V     are  given  bounded  functions  of   e. 

(5) 
Under  the  above  conditions,  Kruskal ^    has  shown  that 

a  formal  solution  of  the  initial-value  problem   (l.l),   (1-2) 

can  be  obtained  by  assuming  for  R(t;€)   a  series  expansion 

of  the  following  kind: 

R^  ^  el"!  R^  e^^  ^/^  (lO) 

n=-aD 

Here  e(t,€)   is  determined  (except  for  an  additive  constant) 
by 

e  =  b[Rq,  ti  €]  (1.4) 


The  coefficients  R   are  obtained  as  power  series  in  £: 

R^(t,^)  ~  ^  Rnic  ^"^  (1-5) 

n=o 

The  vectors  R    may  depend  on  e  as  well  as  on  t  if 

B  and   E  depend  on   e .   To  first  order  in   £,   the  series 


(1.3)   represents  tne  motion  as  a  circular  gyration,  given  by 
the  terms   e  R    e"""^  ,   about  a  moving  guiding  center  given 
by  Rq.   The  frequency  of  gyration  Is   e/e ,   which,  by   (1.4), 
is  the  cyclotron  frequency  evaluated  at  the  guiding  center. 

The  coefficients  R^^  are  obtained  by  Inserting   (1-3) 
Into  the  differential  equations   (l.l)   and  Initial  conditions 
(1.2),  computing  as  If  the  series  did  In  fact  converge,  and 
as  If  B  and  E  were  analytic  functions  of  R.   The  vector 
functions   B,   E  are  expanded  In  Taylor  series  about   R-. 
All  terms  multiplied  by  e"""  '         (for  each  k)   are  collected. 
The  total  coefficient  of  each  exponential   e       Is  set  equal 
to  zero.   This  procedure  yields  equations  for  the  vectors   R  . 
These  equations  are  then  solved  by  assuming  the  expansion   (1-5) 
for  each  R   and  then  solving  the  resulting  equations  --  some 
of  which  are  finite  equations,  and  some  of  which  are  differential 
equations.   Kruskal  ^'^''  has  shown  that  this  procedure  determines 
the  coefficients   R  ,   uniquely,  and  in  sequence.   The  R  i,'s 
for  small   n,   k  are  completely  determined  before  the  R  ,   for 
larger  values  of   n,   k.   In  this  sense,  then,  Kruskal  has  shown 
that   (1-3)   Is  a  formal  solution  of  the  initial-value  problem. 

Our  purpose  in  this  paper  Is  to  show  that   (1-3)   is  in 
fact  an  asymptotic  representation  of  the  exact  solution  of  the 
initial- value  problem,  valid  as   e   tends  to  zero. 
2.   Statement  of  the  Results. 

To  state  our  results  compactly  it  is  convenient  to  use 
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the  notation   "0(€  )"   --   to  be  read  "a  quantity  of  the  order 
of  magnitude  of   e  ."   To  explain  what  this  means,  suppose 
that   f(t;  e)   is  a  quantity  depending  on   t,   e,   and  we  write 

f  =  Oie"^) 

Then  we  mean  that  given   t,  >  0   there  Is  a  positive  number 
A(t-,),   depending  on   t,  ,   but  not  on   €,   such  that 

for   0  <  t  <  t^   ,   |f I  <  A€" 

It  follows  from  Kruskal's  analysis  that  If   N  >  2   Is 

nou) 
(N) 


a  positive  Integer,  and  R^  '   Is  the  sum  of  a  large  enough. 


finite,  number  of  terms  of  the  series   (I.3),   that   R 
approximately  satisfies  the  equations  of  motion  and  Initial 
conditions,  In  the  following  sense: 

6r(N)  +  b(rN,  t;  e)XR^^^  "  E(R^N^  t|  €)  =  o{^^)       (2.1) 


for   t  =  0,   R^^^  =  R  +  O(e^) 
r(N)  =v^0{^'') 


(2.2) 


The  error  estimates  In   (2.1),   (2.2)   may  be  obtained  by 
using  Taylor's  theorem  with  remainder.   We  do  not  wish  to 
repeat  Kruskal's  analysis  here,  so  we  shall  assume  that   (2.1), 
(2.2)   have  been  established.   Assuming  this,  and  assuming  that 
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B,      E,  S,   V.  have  the  properties  stated  In  the  Introduction, 
we  shall  prove  the  following  theorems: 

Theorem  I:   The  Inltlal-value  problem  stated  by   (l.l), 
(1.2)   has  exactly  one  solution,  which  is  defined  for  all 
values  of  t. 

Theorem  II;   This  solution  R,   and  the  approximation 
R^  ',   and  their  derivatives,  are  bounded  uniformly  In  e. 
That  is. 


R  =  0(1)   ,   R  =  V  =  0(1) 


R 


(N)  ^  o(,)   ,   ^(N)  ^  ^(N)  ^  o(,) 


(2.4) 
(2.5) 


M 


Our  main  result  Is  to  the  effect  that   R^  '      not  only 
approximately  satisfies  the  differential  equation   (l-l)   and 
initial  conditions   (1-2)   but  is  indeed  approximately  equal 
to  the  exact  solution  R.   We  state  this  result  as  follows: 
Theorem  III: 


R  -  r(N)  .  o(eN-l) 


V  -  v^N)  ^  o(,N-2) 


(2.6) 
(2.7) 


(when  we  say  that  a  vector  is   0(£'^)   we  mean  that  each  component 
of  the  vector  is   0(e^).)   Theorem  III  is  a  precise  statement  of 
what  is  meant  by  saying  that  Kruskal's  series   (1-3)   is  a  valid 
asymptotic  development  of   R(t;  e)   as   €   tends  to  zero. 


3-   Proof  of  Theorems  I,  II, 

Theorem  I  Is  not  quite  trivial,  and  It  Is  necessary  to 
be  assured  of  Its  truth  to  know  that  It  makes  sense  to  talk 
about  "the  solution"  of  our  Inltlal-value  problem. 

We  write   (l-l)   as  a  first-order  system  as  follows: 

R  =  V  (3.1) 

1   ^v   .      .1 


V  = 


T  ^X^  ^  T^  (5-2) 


The  assumptions  about  the  continuity  and  differentiability  of 
B,   E  guarantee  that,  for  sufficiently  small  values  of   t, 
equations   (3.I),   (3-2)   have  a  unique  solution  satisfying  the 
initial  conditions   (1.2).   Convelvably,  however,  this  solution 
can  fall  to  exist  beyond  a  finite  time.   However,  this  can  happen 
only  if   R_  or  _v   become  unbounded  in  a  finite  time.   (That  is, 
if   |r|   or   |v|   either  tend  to  infinity,  or  oscillate  with  an 
amplitude  which  tends  to  infinity  as   t   tends  to  some  finite 
upper  limit   t„).   Thus,  to  show  that  the  solution  is  defined 
for  all   t,   it  is  sufficient  to  show  that  the  solution  is 
bounded  in  any  finite  time.   We  do  this  by  proving  Theorem  II, 
according  to  which  not  only  are  R,   v   bounded,  but  are  bounded 
by  bounds  which  do  not  depend  on  £. 

To  prove  Theorem  II  we  transform  (3.I),  (3-2)  by 
introducing  in  place  of  v  the  velocity  w  relative  to  a 
coordinate  system  moving  with  the  "drift  velocity"   v^,   given  by 
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D 

Note  that  --  because,  by  assumption,   B~   Is  bounded  uniformly 
In  R,   t,   €   --   Vq  has  all  the  continuity,  differentiability, 
and  boundedness  properties  which  were  assumed  for  B,   E.   We 
define  w,   then,  by 

V  =  Vq  +  w  (3.^) 

Then   (3-2)   becomes 

1         1  (B-E) 
B 

The  derivative  of  v^  can  be  expressed  as  a  function  of 
w,   R,   t   (and  €)   by 

Vq  =  [ ( Vq  +  W)  .  V ]  Vq  +  ^  Vq 

where  'y'   Is  the  gradient  operator  with  respect  to  R.   We  now 
form  the  scalar  product  of  w  and   (3-5),   and  obtain 

d_  (i  w^)  =  -  ^^  (w.B)  -  w.-£^  v„  -  w.([v_+w].Vv^l  (3.6) 
dt   2         e    B^  ~c}t         '-"^      ~^  J 

Here  we  have  used   d/dt.   Instead  of  the  dot,  to  denote  total 
time  differentiation,  and  we  have  used  the  symbol   w   to  denote 
the  length  of  w;   that  is. 


11 


W   =   y/W.W 

Equation   (3-6)   Is  the  equation  for  the  rate  of  change  of 
the  kinetic  energy  of  the  particle,  as  defined  in  the  moving 
coordinate  system. 

The  Important  feature  of   (3-6)   for  our  purposes  is  that 
the  right-hand  side,  considered  as  a  function  of  the  components 
of  w,   is  a  quadratic  polynomial.   The  coefficients  of  this 
polynomial  are  not  known,  but  they  are  known  to  be  bounded. 
This  follows  from  the  assximptions  stated  in  the  introduction. 
(Here  we  make  essential  use  of  the  assumption  that   €   (E*B) 
is  bounded.)   Let   Q(w)   denote  the  right-hand  side  of   (3.6). 
Because   Q(w)   is  a  quadratic  polynomial  with  bounded  coefficients, 
it  follows  that   Q,(w)   is  bounded  as  follows: 

|Q(w)|  <  Aw^  +  c  (3.7) 

where   A,   C  are  independent  of  R,   t,   e.   Hence,  by 
(3.6) ,   we  have 

—  (w^)  <  I  — (w^)|  <  2A  w^  +  2C 

dt         dt 

whence 

^  (e-2At  w2)  <  2Ce-2At 

dt 

so  that,  on  integrating,  we  obtain 
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2  2 

where  w   means  the  initial  value  of  w  .   But   (^-8)   Implies 

w^  =  0(1) 


arid  thus 


w  =  0(1) 


Hence 


V  =  Vq  +  w  =  0(1)  +  0(1)  =  0(1)  0.9) 

and,  since 


R  =  V 


we  have 

R  =  0(1)  (3.10) 

The  same  arguments  can  be  applied  to  the  function 
R^  ^   which  satisfies   (2.1),   (2.2).   The  occurrence  of  an 
additional  term  which  is   0(e    )  (N  >  2)   on  the  right-hand 
sides  of  some  of  the  equations,  does  not  affect  the  results. 
Thus  we  have 
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V 


(N) 

,(N)  _ 


0(1)  (3-11) 


R^^^  =  0(1)  (5-12) 

Having  established   (5-9)  -  (3-12)   we  have  proved  Theorem  II, 
and  at  the  same  time,  have  completed  the  proof  of  Theorem  I. 

4.   Some  Properties  of  the  Vector-Product  Operator. 

To  prove  Theorem  III  we  shall  have  to  make  a  further 
transformation  of  the  equations  of  motion.   We  obtain  this 
transformation  by  decomposing  w_  Into  a  sum  of  eigenvectors 
of  the  vector-product  operator  V  defined  by 

Vw  =  IB  Y  w  (4-1) 

(The  factor  1  In  (4.1)  makes  V  hermltlan) .  Therefore 
we  Insert  here  a  discussion  of  some  of  the  properties  of  V 
which  we  shall  use. 

B^   using   (4.1),   we  can  compute  the  Iterates  of   V,   as 
follows:  / 

v\  =  IB^BX^  =  B^VW 
Thus   V  satisfies  the  characteristic  equation 

V^  =  B^V  (4.2) 
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Thus  the  eigenvalues  of  V  are  0,   B,   -B,   or,  for  short, 
kB  (k  =  -  1,  0,  l) .   Any  vector  y  can  be  written  uniquely  as 
a  sura  of  three  eigenvectors  of  V  associated  with  the  three 
eigenvalues  kB.   That  Is,  for  any  y  there  Is  exactly  one 
set  of  vectors  y,  (k  =  -  1,  0,  l)  such  that 

1 

y  =  Zi  +  Zo  ^  ^-1  "   y       ^k  ^^"^^ 

k+  -  1 


where 


V  y^  =  kB  y^  (4.4) 


The  vector  y.   Is  a  linear  function  of  y,   given  by 


l^  =  \y  (^-5) 


The  projection  operators   P,   are  given  by 

P,  =  At  (V^+BV)  (4.7) 


1    2b2 


P 


1    2b2 


(V^-BV)  (4.8) 


where   I   Is  the  Identity  operator. 

To  prove  these  statements  we  note  the  following  relations 
-  easily  verified  by  direct  computation.  If   (4.2)   Is  used  to 
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express      V   ,      V       in  terms   of     V,      AP : 

V  P^^  =  P^^  V  =  k  B  P^^  (^.9) 

P^P^   =  0     for     k   ^   ^  (^-10) 

Pq  +   Pi   +   P_i    =   I  (^-11) 
By      (4.10),       (4.11)      we  have 

\'  =  Pk  y~  ^i-\^  =  \  (^-12) 


Now   (4.9),   (4.11)   show  that 


k=  -  1 


Is  a  decomposition  of  y   of  the  type  required.   Conversely, 
the  relations   (4.4),   (4.9)   imply 


so  that 

^A  =  °  ^°^   i  +  k  . 
Consequently,  it  follows  from   (4. 3),   (4.11),  that 
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i=  - 1 


so  that   (4.13)   is  the  unique  decomposition  of  y  into 
eigenvectors  of  V. 

3.   Transformation  of  the  Equations  of  Motion. 

We  have  seen  that  the  vector  w  defined  by   (3-4) 
satisfies  the  differential  equation 

w  =  |-  V  w  +  f  (R,  t,  w;  e)  (5-l) 

The  vector  function  f   is,  in  detail,  given  by 

f  =1  (1^)  B-  w.^v^-  (v^+  w).Vv^  (5.2) 

We  note  that   f   Is  a  continuous  function  of   R,   t,   w  with 
continuous  derivatives  with  respect  to   R,   t,   w.   It  follows, 
by  Theorem  II,  that,  for  the  values  of  R,   w  defined  by  the 
solution  of  the  initial- value  problem,   f  and  its  partial 
derivatives  with  respect  to  R,  t,      w  are  all   0(l).   The 
vector  R  satisfies  the  differential  equation 

R  =  Vq(R,  t;  €)  +  w  (5.3) 

The  initial  values  of  R,   w  are  given  by 
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for  t  =  0,    R  =  R 

(5.4) 
w  =  w 

where 


w  =  V  - 


b2 


where   E,   B  are  evaluated  at   t  =  0,  R  =  R. 

The  equations  of  motion  in  the  form   (5-1) j   (5-3) 
are  not  well  adapted  to  the  study  of  the  behaviour  of  the 
solution  as   e   tends  to  zero.   The  reason  lies  in  the  factor 
/€   occurring  in  the  right-hand  side  of   (5-1) •   A  Lipschitz 
constant  for   (5-1)   would  be  of  the  order  of   /e .   Proceeding 
directly  from   (5-1)   to  estimate   R.  -  R     one  would  obtain  an 
estimate  which  would  tend  to  infinity  as   e   tended  to  zero,  and 
which  would  hence  be  useless. 

We  avoid  this  difficulty  by  a  transformation  which  begins 
by  decomposing  w   into  eigenvectors  of  V.   We  have 


1 

w  =  ^     w^  (5-5) 

k=  -  1 


where 


— k  "  ^k  -  (5-6) 


Vw^  =  kBw^  (5-7) 
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Differentiation  of   (5-6)   shows  that 

^k  ^  =  ^k  -  ^k  ^  (5.8) 


Thus  If  one  operates  on   (5-1)   with  P,   one  obtains 

w,  -^«,  =  0^  (5.9) 

Where 


^^  =  \^+  \L 


We  write  G,   as  a  function  of  w,  ,   R,   t,   e   by  writing 


Pk  =  ( Vq  +  w)  .  V  P^  +  ^  P^ 


1 


and  then  inserting   \      w,   for  w   everywhere. 

k=  -  1 


The  next  step  is  to  introduce  new  dependent  variables 
u,   by  the  definition 

u^  =  e-^^^/^  w^  (5.10) 

where  the  function   §(t)   is  defined  by 


t 
i(t)  =  /  B(R(t',  e),    t';  G)dt'  (5.11) 

0 
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Then   (5-9)   becomes 

K   =  ^      ^         9.1^    (LLRj^^^         ,  R,  t;  e)  (5.12) 

Now  we  make  a  further  transformation.   We  introduce  a  new 
Independent  variable  )6     and  a  new  dependent  variable  T{j6). 
The  function  T{)6)      Is  defined  as  the  Inverse  function  of  the 
function   J(t)   defined  by   (5.11).   (This  Inverse  function 
exists  and  is  monotone  increasing  since  because   B  >  Bnln,   the 
function  I      in   (5-11)   is  monotone  increasing.)   That  is, 
T{^)      is  defined  by  the  relations 

$  {T{0))    =  0 

(5.13) 
T  (I(t))  =  t 

Thus  in   (5.12)   we  can  replace  $     tiy     0     and   t   by   T.   ESy 
(5.11)   we  have 

dT  _  1 (^   ih\ 

00  ~   B(R(T:  e),  T,  e)  ^^    '' 

o 

By      (5.12)   we  have,  consequently. 


S^^ie-^^^/^O,  (tu,e^^^/SH,  T,  .)  (5.15) 


and   (5.3)   becomes 
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t   =  i   v^(H>  T;  e)  +  1  V-  u^e^^^A  (5.16) 


Let  X  be  a  vector  with  a  large  number  of  components  - 
(thirteen)  -  the  components  of   X  being  the  scalar   T,   the 
components  of  R,   and  all  the  components  of  the  vectors   u,  , 
arranged  in  some  definite  order.   Then  according  to   (5-14), 
(5- 16),   (5-15),   the  vector  X  satisfies  a  differential 
equation  of  the  following  form 

^  X  =  P  (X,  jz^;  €)  (5.17) 

Note  that  the  vector  function  P   is  a  continuous  function  of 
the  components  of   X,   having  continuous  partial  derivatives 
with  respect  to  the  components  of   X.   Purthermore,   F  and  its 
partial  derivatives  are  bounded  independently  of   e   in  the 
following  sense:  if   t,  >  0   there  is  a  number  K(t|)   such  that  if 

\0\    <  Bmin  t,  (5.18) 

then 

|P|  <  K 

(5.19) 
<  K 

where  P'   means  any  derivative  of  P  with  respect  to  a  component 
of  X,   and  the  notation   |p|   means  the  sum  of  the  absolute  values 
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of  the  components  of  F.   To  see  this,. note  that   (5-1^)   and 
(5.18)   imply 

|T|  <  t, 

so  that   (5.19)   follows  from  Theorem  II  and  the  assumptions 
listed  in  the  Introduction. 

The  vector  X  satisfies  initial  conditions,  which  can 
be  inferred  from   (5-^),   to  the  effect  that 

for    JZ^  =  0,  X  =  X  (5-20) 

(including  the  condition 
for    0=0,    T  =  0. ) 

A  similar  transformation  can  be  performed  on  the 
equations   (2.1),   (2.2)   which  govern  the  putative  approximate 
solution  R^  ' .   We  define 

i(^^(t)  =  /   B(R(^^(t',  e);  t'€)dt'  (5.21) 

0 

(N)  _   -iki^N)/,   (j^)  (5.22) 

-iik    ~  -k 

^N)  ^    ^(N)  (5.23) 

— k      k  — 

(Here   P,   is  regarded  as  a  function  of   B,   which  in   (5-23)  is 
a  function  of  R^^' ,  t,  e.)   Also,  we  define  the  function  ^ 
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(5.24) 

t(n)  (/(N)(t))  =  t 

We  find  that  the  corresponding  vector  X^^   satisfies  the 
differential  equation 

^^W    .p.  (x(N),  0;    e)   +f(eN-l)  (5.25) 

and  the  Initial  conditions 
For    0=0,      X^^^    =  X  +  O(eN)  (5-26) 

The  differential  equations   (5. 17),   (5-25)   with  the 
Initial  conditions   (5.20),   (5-26)   will  be  used  to  complete 
the  proof  of  Theorem  III. 

6.   Proof  of  Theorem  III. 

It  Is  clear  that  by  the  same  reasons  which  were  given 
for   (5.19)   that   K(t,)   can  be  found  such  that 

for    \0\    <  Bmin  t,  (o-l) 

we  have 

IZi  <  K 
|F' I  <  K 
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if  F,   F'   are  evaluated  anywhere  on  the  line  segment  joining 
X{0)      and  X   (-(2^)  •   It  follows,  then,  that 

|F(X,  0;   ^)    -  F   {X^^K  0;   ^)\    <K   \X-   X^^^  I      (6.2) 


If   (6.1)   is  satisfied.   The  inequality   (6.2)   states  that 
F  satisfies  a  Lipschltz  condition  uniformly  in  € .   It  is 
now  an  easy  matter  to  compare  the  vectors  X{0) ,    X   (0) . 

To  do  this,  we  subtract  the  differential  equation 
(5-25)   for  X^  from  the  equation   (5-17)   for  X;      using 
(6.2)   we  obtain 

1^  (X  -  X^^hl    <  K  |X  -  X^^^l  +  0  (€^-1)         (6.3) 

If  we  subtract  the  initial  conditions   (5-26)   from  the 
initial  conditions   (5.20),   we  obtain 

for  0=0,      |X  -  X^^^l  =  0  (e^)  (6.4) 


Let  us   now   integrate      (6.3);      we   obtain 

|X-    X^l    -    0    (cN)    <    I    /^    (X-    X^^^ldj^l 


(6.5) 


<  /  \m  (^  -  x^^^)l  djz^  <  K  /|x  -  x^^^\a0^  o(€N-i) 

0        ^-"^  0 


If      H      is   defined  by 
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H  =   /   |X  -  X^^^l  d0 
0 

then   (6.5)   implies 

^  <  KH  +  O(eN-l)  (g_g^ 

Noting  that 

for   0=0,  H  =  0 

we  conclude  from   (6.6),   as  in  a  similar  calculation  we 
have  done  before,  that 


H  =  0(£N'^) 


and  hence,  by   (6.5),   we  have 

X  -  /  =  0(e^"^)  (6.7) 

Now  in  (6.7),  the  vectors  X,  X^  ''  are  compared  at 
the  same  value  of  ^.  What  we  must  do  to  obtain  Theorem  III 
is  to  compare  R,  R^  '  and  u  ,  u  at  the  same  values 
of   t^.   We  now  procede  to  make  this  comparison. 

We  have,  if   0  =  l^^' (t), 

i  > 

-  |R(t)  '-  R(N^(t)|  =  |R(T^^^0))  -  r(^^(t(^^(0))|    (0.8) 
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and 

|R(t(^\0))  -  r(N)(t(N)(^))|  <  |R(T(iZi))  -  R(NhT^^^0))| 

+  |R(T(0))  -  R(T(N^iZS))|  (6.9) 

Now   (6,7)   implies.  Inter  alia,  that 

|R(T(0))  -  r(N)(t(N)(0))|  =  o(e^-^)  (6. 10) 

and 

|T(0)  -  t(N)(0)|  =  0(€N-1)  ^g^^^j 

By      (5-3)   we  see  that 

R  =  0(1) 
and  hence,  using   (6. 11),   we  see  that 

|R(T(0))  -  R  (t(N)(0))|  <  0(1)  .  0(€N-1) 
so  that 

R  (T(0)  -  R  (T^^^0))  =  0(6^"^)  (6.12) 

Combining   (6.8),   (6.9),   (6. 10),   (6.12)   we  conclude  that 

R  -  R^N^  =  0(eN"^)  (6.13) 
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Likewise,  noting  by   (5.12)   that 

u^  =  0(1) 
we  obtain 

By   (6.15)   and  the  definitions   (5-11),   (5-21)   of   f,  i^^', 
we  see  that 

Hence,  by  the  definitions  of   u^,   u^^  ,   we  obtain  (using 
Theorem  II) 

I   N,   I  iki/e   ^ik^^^V^I 
+  |u^  I   |e   /   -  e      '     \ 

<  0  (€^-1)  +  0(1)  .  I  .  0(€N-1)  ^  o(,N-2)  (5_15) 

Adding  the  inequalities   (6.15),   we  have,  since 

(N)  _T—  .  (N) 


w  =  ■.    w^,  w^  '  =  \   w, 


the  result 
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w  -  wW  =  o(6N-2) 


whence,  using   (6.15),   we  see  that 

V  -  vW  =  o(€N-2)  (g  ^^) 

Now  that   (6.13)   and   (6.16)   have  been  established. 
Theorem  III  Is  proved,  and  our  task  Is  finished. 
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